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OO ! Abstract 

Vbfnlo is a fully flexible parton level Monte Carlo program for the simulation of 
vector boson fusion, double and triple vector boson production in hadronic collisions 
^ ■ at next-to-leading order in the strong coupling constant. Vbfnlo includes Higgs and 

vector boson decays with full spin correlations and all off-shell effects. In addition, 
Vbfnlo implements CP-even and CP-odd Higgs boson via gluon fusion, associated 
with two jets, at the leading-order one-loop level with the full top- and bottom-quark 
mass dependence in a generic two-Higgs-doublet model. 

A variety of effects arising from beyond the Standard Model physics are im- 
plemented for selected processes. This includes anomalous couplings of Higgs and 
vector bosons and a Warped Higgsless extra dimension model. The program offers 
the possibility to generate Les Houches Accord event files for all processes available 
at leading order. 
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Program summary 



Program Title: Vbfnlo 
Journal Reference- 
Catalogue identifier: 
Licensing provisions: GPL version 2 



Program obtainable from: http : //www-itp . particle . uni-karlsruhe . de/vbf nlo/ 



Distributed format: tar gzip file 

Programming language: Fortran, parts in C++ 

Computer: All 

Operating system: Linux, should also work on other systems 

Keywords: NLO Monte Carlo program, one-loop QCD corrections, electroweak bosons, hadronic 
collisions 

PACS: 11.15.-q, ll.80.Cr, 12.38.Bx, 12.60.Pr 
Classification: 11.1, 11.2 

External routines/libraries: Optionally Les Houches Accord PDF Interface library and the GNU 
Scientific library. 

Nature of problem: To resolve the large scale dependence inherent in leading order calculations 
and to quantify the cross section error induced by uncertainties in the determination of parton 
distribution functions, it is necessary to include NLO corrections. Moreover, whenever stringent 
cuts are required on decay products and/or identified jets the question arises whether the scale 
dependence and a k-factor, defined as the ratio of NLO to LO cross section, determined for the 
inclusive production cross sections are valid for the search region one is interested in. 
Solution method: The problem is best addressed by implementing the one-loop QCD corrections 
in a fully flexible NLO parton-level Monte Carlo program, where arbitrary cuts can be specified 
as well as various scale choices. In addition, any currently available parton distribution function 
set can be used through the LHAPDF library. 

Running time: Depending on the process studied. Usually from minutes to hours. 
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1 Introduction 

The physics potential of the TeVatron and even more of the Lhc relies, to a large 
extent, on our ability to provide accurate cross section predictions both for signal and 
background processes. The latter are often generated by QCD interactions followed by 
weak transitions of partons to vector bosons. A precise description of such hard QCD 
production processes is needed, as well as a method for simulating the measurable hadronic 
final states. Reaching these goals requires next-to-leading order (NLO) QCD calculations 
presented in the form of parton level Monte Carlo (MC) generators which are an efficient 
solution when it comes to final states characterized by a high number of jets and/or 
identified particles. When kinematical cuts are imposed, as is mandatory for processes 
involving QCD radiation, analytical phase space integration becomes impractical and 
implementation of results in the form of Monte Carlo programs becomes the method of 
choice. 

Vbfnlo is a fully flexible MC program for vector boson fusion (VBF), double and 
triple vector boson production processes at NLO QCD accuracy. Since real emission 
processes are part of the NLO cross sections, Vbfnlo provides the possibility to calculate 
cross sections for the corresponding process with one additional jet at leading order (LO) 
in the strong QCD coupling. In addition, the simulation of CP-even and CP-odd Higgs 
boson production in gluon fusion, associated with two additional jets, is implemented at 
leading order in the strong coupling with the full top- and bottom-quark mass dependence 
in a generic two-Higgs-doublet model. Several models for anomalous couplings of Higgs- 
and vector bosons and a Warped Higgsless extra dimension model have been implemented. 

Arbitrary cuts can be specified as well as various scale choices. Any currently available 
parton distribution function (PDF) set can be used through the LHAPDF librarjill. For 
processes implemented at leading order, the program is capable of generating event files 
in the Les Houches Accord (LHA) format [20]. 

^ http : //pro j ects . hepf orge . org/lhapdf /| 
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2 Processes 



In the following sections, we describe all production processes and decay modes imple- 
mented in Vbfnlo, together with references to a more detailed discussion of the under- 
lying calculations. 

In the phase space regions which are accessible at hadron colliders, VBF reactions 
are dominated by t-channel electroweak gauge boson exchange. In Vbfnlo therefore 
s-channel exchange contributions and kinematically suppressed fermion interference con- 
tributions [1-3] are disregarded. Throughout, we consider only decays of the weak bosons 
into different lepton generations, such as ZZ i~ i'~ . Results for leptonic final 

states with any combination of leptons (e.g., i~^i~i~^£~) can be obtained thereof by mul- 
tiplying the respective results with appropriate combinatorial factors. Numerically small 
contributions from Pauli-interference effects for identical charged leptons are disregarded. 

2.1 VBF Higgs production in association with two jets 

Hjj production via VBF mainly proceeds via electroweak quark-quark scattering pro- 
cesses like qq' — > qq'H and crossing related reactions. In Vbfnlo, tree level cross sections 
and NLO QCD corrections to the t-channel production process are provided. The subse- 
quent decay of the Higgs boson is simulated in the narrow width approximation (NWA) . 
For the H — > W~^W~ and the H ZZ modes, full off-shell effects and spin correlations 
of the decay leptons are included. The available production process and the decay modes 
are listed with the corresponding process IDs in Table [U Details of the calculation can 
be found in Ref. [4]. 



ProcId 


Process 


100 


PP - 


-Hjj 






101 


PP - 


-Hjj- 


11 jj 




102 


PP - 


-Hjj- 


^"^/^~ jj 




103 


PP - 


-Hjj- 


^ r+T~ jj 




104 


PP - 


-Hjj- 


bbjj 




105 


PP - 


-Hjj- 


w+w- jj 




106 


PP - 


-Hjj- 


^zzjj^e+t 


-£'+£'- jj 


107 


PP - 


-Hjj- 


^ ZZjj ^ i+t 


'Vi'Uli jj 



Table 1: Process IDs for pp Hjj production via weak boson fusion at NLO QCD 
accuracy. 



2.2 VBF Higgs production in association with three jets 

Adding an extra parton to the Higgs production processes of Sec. 12.11 gives rise to Hjjj 
final states. The corresponding cross sections are implemented at NLO QCD accuracy in 
Vbfnlo. A list of all available modes and corresponding process IDs is given in Table [2l 
Details of the calculation can be found in Ref. [5]. 
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1 1 fl 


pp - 


H jjj 






111 

± ± ± 


pp - 


H jjj - 


11 jjj 






pp - 


H jjj - 


-* Z^'*'^" jjj 




1 1 Q 


pp - 


H jjj - 


^ r+T" jjj 






pp - 


H jjj - 


-> bbjjj 




111^ 
1 ID 


pp - 


H jjj - 


^W+W-jjj- 


i-^vil' v^, jjj 


116 


pp - 


H jjj - 


-^ZZjjj^i+ 


i-i'+i'-jjj 


1-1-1 


pp - 


H jjj - 


-^ZZjjj^i+ 


^'I'e'^e'jjj 



Table 2: Process IDs for pp —>■ Hjjj production via weak boson fusion at NLO QCD 
accuracy. 



ProcId 


Process 


120 


pp - 


Zjj 


i+rjj 


121 


pp - 


Zjj 


ViVi jj 


130 


pp- 


^W+jj 


^"^vijj 


140 


pp- 


^W-jj 


^~vaj 



Table 3: Process IDs for Zjj and W jj production via weak boson fusion at NLO QCD 
accuracy. 



2.3 VBF production of a vector boson and two jets 

Vector boson fusion processes can also produce final states with two leptons plus two jets, 
which are generically referred to as "VBF Zjj and W^jj production". These reactions are 
implemented to NLO QCD accuracy in Vbfnlo, see Table [3l Details of the calculation 
can be found in Ref. [6]. 

2.4 VBF production of two vector bosons and two jets 

The production of four leptons plus two jets in the final states at order 0{a^) is dominated 
by VBF contributions. In Vbfnlo, all resonant and non-resonant t-channel exchange 
contributions giving rise to a specific leptonic final state are considered. For simplicity, 
we refer to these reactions as "VBF W~^W~jj, ZZjj, and W^Zjj production". Finite 
width effects of the weak bosons and spin correlations of the decay leptons are fully 
retained. The available processes and corresponding process IDs are listed in Table [H 
Details of the calculation can be found in Refs. [7-9]. 
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T-* T? C C 


zuu 


- 


w+w-jj 






- 




V P- Pl+ nn 
I I I JJ 




pp - 


ZZ jj t 




220 


pp - 


^ W+Zjj ^ 


l+u^l'+l'- jj 




pp - 


-> W-Zjj ^ 


i-ue£'+f'~jj 



Table 4: Process IDs for W~^W jj , ZZjj and W^Zjj production via weak boson fusion 
at NLO QCD accuracy. 



ProcId 


Process 


300 


pp - 


w+w- 


— > i 




400 


PP - 


w+w- 


z - 




410 


PP - 


ZZW+ 






420 


pp - 


zzw- 






430 


pp - 


w+w- 


w+ 




440 


PP - 




w- 





Table 5: Process IDs for the W^W , WWZ , ZZW and WWW production processes at 
NLO QCD accuracy. 



2.5 Double and triple vector boson production 

The production of four- and six-lepton final states mainly proceeds via double and triple 
vector boson production with subsequent decays. In Vbfnlo, the processes listed in 
Table [5] are implemented to NLO QCD accuracy, including full off-shell effects and spin 
correlations of the final state leptons. Details of the calculation can be found in Refs. [10, 
11]. 

2.6 Higgs production in gluon fusion with two jets 

CP-even and CP-odd Higgs boson production in gluon fusion, associated with two addi- 
tional jets, is a process which first appears at the 1-loop level which, therefore, is counted 
as leading order in the strong coupling. This process is simulated including the full mass 
dependence of the top and bottom quark running in the loop of a generic two-Higgs- 
doublet model. The relevant process ID is given in Table [6l Details of the calculation can 
be found in Refs. [12-16]. 
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JT KUL'IIjC'C) 


41 nn 


pp ^ n jj 



Table 6: Process ID for the LO Higgs plus 2 jets production via gluon fusion. 



3 Installing VBFNLO 

The source code of the current version of Vbfnlo can be downloaded from the Vbfnlo 
web-page 

http : //www- itp .particle .uni-karlsruhe . de/vbfnlo/ 

and includes a GNU conforming build system for portability and an easy build and in- 
stallation procedure. 

3.1 Prerequisites 

The basic installation requires GNU make, a Fortran77@ and a C++ compiler. Vbfnlo 
offers the possibility to use the LHAPDF0 library for parton distribution functions. In 
case the simulation of Kaluza-Klein resonances should be enabled, an installation of the 
GNU Scientific Library (GSL)0 is required. 

3.2 Build and installation 

After unpacking the source archive and entering the source directory, the configure script 
can be invoked with several options, a complete list being available via ./configure 
-help. Among these, the most important ones are: 

• --pref ix= [path] 

Install Vbfnlo in the location given by [path] . 

• --enable-processes= [list] 

By default, the code for all available processes is compiled. Optionally, [list] gives 
a comma-separated list of selected processes to be compiled. Possible process names 
are: 

vbf Vector boson fusion processes 

diboson Double gauge boson production 
triboson Triple gauge boson production 

hj j j Higgs boson plus three jet production in vector boson fusion 

ggf Higgs boson plus two jets via gluon fusion 

^The following compilers have been tested: g77 and gf ortran 
" http : //proj ects .hepf orge . org/lhapdf / 
' http : //www . gnu . org/ sof tware/gsl/| 
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• --disable-NLO 

Disable the next-to-leading order QCD corrections. With this option, compilation 
time is shortened. 

• --enable-kk 

Enable simulation of Kaluza-Klein resonances. Disabled by default, the Kaluza- 
Klein option requires the installation of the GNU Scientific Library to be specified 
via --with-gsl. 

• --with-LHAPDF=[path] 

Enable the usage of LHAPDF instead of the built-in PDF sets. Disabled by default, 
[path] specifies the location of the LHAPDF installation. 

• --with-gsl= [path] 

Enable usage of the GNU Scientific Library, [path] specifies the location of the 
GSL installation. 

Once configure finished successfully, the make and make install commands will 
compile and install Vbfnlo, respectively. 

3.3 Special notes for Mac OS X 

Due to a linker problem on this platform, Vbfnlo has to be compiled with static libraries 
disabled by adding the --disable-static flag to the call of the conflgure script. For the 
dynamic libraries to be properly resolved, the environment variable DYLD_LIBRARY_PATH 
has to be set to the library directory of the installation: 

export DYLD_LIBRARY_PATH= [prefix] /lib/VBFNLO 

where [prefix] is the installation directory as chosen by the --prefix parameter. 

3.4 Source and installation directory layout 

The Vbfnlo source tree contains the following subdirectories: 

• amplitudes/ : Routines to calculate matrix elements for the processes provided. 

• doc/ : The source of this manual. 

• helas/ : Helas [17] subroutines used to calculate helicity amplitudes. 

• loops/: One-loop tensor integrals up to flve-point functions. 

• PDFsets/: Built-in parton distributions (CTEQ6L1 and CTEQ6M, [18]). 

• phasespace/ : Specialized phasespace generators for the processes provided. 

• src/ : Source code of the main programs and input files. 

• utilities/ : Routines for administrative tasks, cuts, scale choices and interfaces. 
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The source does not need to be modified to change the simulation parameters. Vbfnlo 
off'ers several kinematical cuts and scale choices. This is illustrated in Sec. [H In addi- 
tion, it provides a few basic histograms. Cuts, histograms and scale choices not al- 
ready provided may be added in the utilities/cuts . F, utilities/histograms . F and 
utilities/scales . F files. 

The installation is performed in a standard UNIX-layout, i.e. the directory specified 
with the --prefix option of the configure script contains the following directories: 

• bin/ : vbfnlo and ggf lo executables. 

• include/VBFNLO/ : Vbfnlo header files. 

• lib/VBFNLO/ : Vbfnlo modules as dynamically loadable libraries. These can also 
be used independently from one of the main programs. 

• share/VBFNLO/ : Input files and internal PDF tables. 

3.5 Running the program 

Both the vbfnlo and ggflo executables contained in the bin directory of the installation 
path do look for input files in their current working directory. An alternative path to input 
files may be specified explicitly by passing the input=path argument to the programs, 
with path denoting the full path where input files are located. 

The input files contained in the share/VBFNLO directory are meant to represent default 
settings and should not be changed. We therefore recommend to symbolically link the 
desired executable and copy the input files to a separate directory. Here, special settings 
may be chosen in the input files and the program can be run in that directory without 
specifying further options. 

3.6 Bug reports 

Please report any problems to 

vbf nloQparticle .uni-karlsruhe . de 

including a short report with which configure options Vbfnlo has been built, as well as 
the versions of compilers and external libraries used. 

3.7 License 

Vbfnlo is distributed under the GNU General Public License (GPL) version 2. This 
ensures that the source code will be available to users, grants them the freedom to use 
and modify the program and sets out the conditions under which it can be redistributed. 
However, it was developed as part of an academic research project and is the result of 
many years of work by the authors, which raises various issues that are not covered by the 
legal framework of the GPL. It is therefore distributed together with a set of guideline^, 
which originally have been formulated and agreed on by the MCnet collaboration for event 
generator software. 

^These guidelines are contained in the GUIDELINES file distributed with the release. 
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4 Input files and parameters 



Vbfnlo is steered through the following input files: 

• vbfnlo.dat: General parameters for a run. 

• ggflo.dat: Additional parameters for the ggflo program. 

• cuts.dat: Values for kinematical cuts. 

• anom_HVV.dat: Parameters for anomalous Higgs couplings. 

• anoiii_WW.dat: Parameters for anomalous triple gauge boson couplings. 

• random . dat : Seed for the random number generator. 

The following subsections will give a detailed description of all available parameters. 

4.1 vbfnlo.dat — general parameters 

• PROCESS : Process ID as described in Sec. H 

• LOPROCESS_PLUS_JET: If set to true, the leading order process with one additional 
jet is generated, i.e. only the real radiation contribution is generated. This option 
is available for all but gluon fusion processes. 

• LEPTONS : Choice of the final state leptons according to the MC particle numbering 
scheme [19]. If the selected configuration is not available, default values are used. 

• LO_ITERATIONS : Sets the number of iterations for the integration of LO cross sec- 
tions. Usually more than one iteration is used in order to adapt the integration grid 
and thus improve the efficiency of the MC integration algorithnj^. For an adapted 
grid file (see LO_GRID) this parameter can be set to 1. Default is 4. 

• NLO.ITERATIONS: Analogous to LO.ITERATIONS, but for the real emission part of 
an NLO calculation. Since the corresponding phase space is different from the LO 
configuration, a second independent MC integration has to be performed. Default 
is 4. 

• LO_POINTS: Determines the number of phase space points that are generated in 
each iteration. In the last iteration there are 2^ points, where A^=LO_POIWTS. In 
each previous iteration, the number of points is half the value of the following one. 
Example: For 4 iterations (LO.ITERATIONS = 4) and LO.POINTS = 20, there are 
2^^ generated points in the first, 2^^ in the second , 2^^ in the third and 2^° ^ 10^ 
in the last iterationQ. Default is = 20. 

• NLO.POINTS: Similar to LO.POIWTS, but for the real emission part of a NLO calcu- 
lation. 

^For all NLO calculations the virtual contributions are calculated using the already optimized leading 
order grid. 

^The virtual contributions are calculated for 2^ points only. 
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• LO_GRID : Sets the name of the grid files that are generated at the end of each itera- 
tion. Choosing name as input parameter, in each iteration X a grid file name . out .X 
will be produced. If a grid file name is already present in the working directory, the 
program reads in this file when executed. 

• NLO_GRID: Similar to LO_GRID, but for the real emission part of a NLO calculation. 

• NLO_SWITCH: Switch for the NLO part of a process, if available. If set to .true., 
cross sections and histograms are calculated to NLO QCD accuracy. Default is set 
to .false.. 

• ECM: The center-of-mass energy y/s of the collider, measured in GeV. Default is 
14000 GeV. 

• BEAMl, BEAM2: Define the type of particle of each beam. Possible options are +1 
for proton beams and —1 for anti-proton beams. Default is proton-proton collisions, 
(+1, +1). 

• ID_MUF: Choice of the factorization scale. See Table [7] for a list of available options. 
Default is 0. 



• ID_MUR: Choice of the renormalization scale. See Table E] for a list of available 
options. Default is 0. 



ID.MUF 


Factorization Scale 


Process class 





user defined constant scale set by MUF_USER 


all 


1 


momentum transfer of exchanged W/Z boson 


vbf 


2 


min(pT(ji),PT(j2)) 


vbf 


3 


invariant mass of VV system 


diboson 


4 


invariant mass of VVV system 


triboson 


5 


\/Pt(Ji) X Pt(J2) 





Table 7: Factorization scale options. 



• MUF_USER: If 1D_MUF is set to 0, this parameter sets the user defined constant 
factorization scale measured in GeV. Default is 100 GeV. 

• MUR_USER: If 1D_MUR is set to 0, this parameter sets the user defined constant 
renormalization scale measured in GeV. Default is 100 GeV. 

• XIF : Factor by which the factorization scale is multiplied. May be used to analyze 
the scale dependence of differential cross sections. Default is 1.0. 

• XIR: Factor by which the renormalization scale is multiplied. May be used to 
analyze the scale dependence of differential cross sections. Default is 1.0. 
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ID MUR 


R P.NORM A T J7 ATTON SCAT.P, 


X JTVi Vj* 12j kJ kJ 'w'J-^jT.OO 





Uoci LiuiiiicLi cuiio Ldiii oL-dic L uy riun_uoEin 


all 


1 


inompntiim tr^msfpr of P'vrVi^iTiP'prl lioson 

X J.AV_/ AAA v^±A D l_l J.A1 D A <Jj Al kJ A A KJL Kj./\.Kj±±(AiL1 S-^KjKA. r r / i-/ KJKJiJKJLl 


vbf 


9 


Ililll l^/^T Ul J 5 FT U2 J J 


vhf 

V UJ. 


Q 
O 


lllVdlldjlll llldoo Ui V V byolclil 


QlDUbUll 


4 


invariant mass of VVV system 


triboson 


5 




6&-^ 



Table 8: Renormalization scale options. 



4.2 vbfnlo.dat — physics parameters 

• HMASS: Standard Model Higgs boson mass in GeV. Default value is 120 GeV. 

• TOPMASS: Top quark mass in GeV. Default value is 172.4 GeV. 

• BOTTOMMASS : Bottom quark pole mass in GeV, used in the calculation of the Higgs 
width and branching ratios as well as in the heavy quark loop in the gluon fusion 
process. Default value is 4.855 GeV. 

• CHARMMASS : Charm quark pole mass in GeV used in the calculation of the Higgs 
width and branching ratios. Default value is 1.65 GeV. 

• ALFA_S : Strong coupling constant as used in the calculation of W, Z and H widths. 
The strong coupling constant used in the matrix element calculations is printed out 
during run time. Default value is 0.1176. 

• EWSCHEME: Sets the scheme for the calculation of electroweak parameters. A sum- 
mary of the four available options is given in Table [9l Default value is 3. 

• FERMI_CONST: Fermi constant, used as input for the calculation of electroweak pa- 
rameters. Default value is 1.16637 x 10~^ GeV~^. 

• ALFA: Fine structure constant, used as input for EWSCHEME = 1 and 4. Within the 
other schemes this parameter is calculated. Default value is 7.2973525376 x 10~^. 

• SIN2W: Sinus squared of the weak mixing angle. Used as input for EWSCHEME = 2 
and 4. Within the other schemes this parameter is calculated. Default value is 
0.23119. 

• WMASS: W boson mass in GeV. Default value is 80.398 GeV. 

• ZMASS: Z boson mass in GeV. Default value is 91.1876 GeV. 

• ANOM_CPL: Options for anomalous Higgs and gauge boson couplings. These are 
available for the Hj] and W~^W~jj production processes in VBF. Default is set to 
.false . . 
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EWSCHEME 

J_j VV kJ v./! lJ_ji iJ_j 


Par amp.tp.r 


r)P,PATIT,T VATJIP 


TnPTIT / G a PPT it, ATP; P) 


1 


FERMI_CONST 
ALFA 
SIN2W 
WMASS 
ZMASS 

£-11 iri. !_/ t~J 


1.16637 X 10-^ GeV"^ 
7.2973525376 x 10^^ 
0.23110323 
79.9595 GeV 
91 1876 GeV 


Input 

Input 
Calculated 
Calculated 

Input 


2 


FERMI.COWST 
ALFA 
SIM2W 
WMASS 
ZMASS 


1.16637 X 10~^ GeV~2 
7.7602239787 x 10"=^ 
0.23119 
79.9544 GeV 
91 1876 GeV 


Input 
Calculated 

Input 
Calculated 

Input 


3 


FERMI.COMST 
ALFA 
SIM2W 
WMASS 


1.16637 X 10-5 GeV-2 
7.5562544251 x 10'^ 
0.22264585 
80.3980 GeV 
91 1876 GeV 

X . X w 1 \J V.' v_/ V 


Input 
Calculated 
Calculated 

Input 

Input 


4 


FERMI.COMST 
ALFA 
SIM2W 
WMASS 
ZMASS 


1.16637 X 10"5 GeV-2 
7.2973525376 x 10"=^ 
0.23119 
80.398 GeV 
91 1876 GeV 


Input 
Input 
Input 
Input 
Input 



Table 9: Electroweak input parameter schemes. 



• KK_MOD : Option for the Warped Higgsless Model. It is available for all VVjj pro- 
duction modes in VBF. Default is set to .false.. 

4.3 vbfnlo.dat — parameters for event output 

Vbfnlo generates parton level events according to the most recent Les Houches Accord 
(LHA) format [20] for processes available at leading order. 

• LHA_SWITCH: Switch on or off output of LHA event files. Default is set to .false.. 

• UMWEIGHTIMG_SWITCH: Option for event weights. If set to .true., events are un- 
weighted (event weight = +1). If set to .false., events are weighted. Default is 
set to .false . . 
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• PREMEVUMW: The number of events used in the last iteration in order to calcu- 
late/estimate the premaximal weight which is needed in the first step of the un- 
weighting procedure. Default is 1000. After all events are unweighted, the maximal 
weight is again calculated and a reweighting procedure is applied. 

• TAUMASS : Option to include the mass of r leptons in the LHA file. So far this option 
only works for the vbf processes. Default is set to .false.. 

4.4 vbfnlo.dat — PDF parameters 

Vbfnlo may use built-in parton distribution functions (PDF) or the LHAPDF library. 

• PDF.SWITCH: Option to choose which PDFs are used. If set to 0, built-in PDFs 
(CTEQ6L1 for LO and CTEQ6M for NLO calculations) are used [18]. If set to 1, 
an interface to LHAPDF is provided via LHAGLUE [21]. 

• LO.PDFSET: LHAGLUE number for the LO PDF set, see PDFsets. index or Ref. [21]. 
Default is 10042 (CTEQ6L1). 

• WLO.PDFSET: LHAGLUE number for the NLO PDF set, see PDFsets. index or 
Ref. [21]. Default is 10000 (CTEQ6M). 

4.5 vbfnlo.dat — parameters for histograms 

Vbfnlo provides output for histograms in the following formats: Topdrawef|§, Root@ 
and GNUPLOT0. 

• ROOT: Enable output of histograms in RoOT format. Default is set to .false.. 

• TOP: Enable output of histograms in TOPDRAWER format. Default is set to .false.. 

• GNU: Enable output of histograms in Gnuplot format. Default is set to .true.. 

• REPLACE: Switch to overwrite existing histogram output files. Default is set to 
.true . . 

• ROOTFILE: Name of the RoOT output file. Default is histograms. 

• TOPFILE: Name of the TOPDRAWER output file. Default is histograms. 

• GNUFILE: Name of the Gnuplot output file. Default is histograms. 

4.6 cuts.dat — parameters for kinematical cuts 

Jet-specific cuts: 



• RJJ_MIN: Minimum separation of two identified jets, ARjj = ^ Ayjj + A0^^-, used 
by the k± jet finding algorithm [22] which combines all partons. Default is 0.8. 




^ http : / / www . pa . msu . edu/ ref erence/topdrawer-docs/ 



http :/ /root . cern . ch/ 



ic 



http : //www. gnuplot . info/ 
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• Y_P_MAX: Maximum allowed pseudorapidity for final state partons. Default is 5.0. 

• PT_JET_MIN: Minimum transverse momentum for identified jets. Default is 20 GeV. 

• Y_JET_MAX: Maximum allowed rapidity for identified jets. Default is 4.5. 

Lepton specific cuts: 

• Y_L_MAX: Maximum pseudorapidity for charged leptons. Default is 2.5. 

• PT_L_MIN : Minimum transverse momentum for charged leptons. Default is 10 GeV. 

• MLL_MIN : Minimum invariant mass for any combination of oppositely charged lep- 
tons. Default is 15 GeV. 

• RLL_MIM: Minimum separation of charged lepton pairs, ARu. Default is 0. 

• RLL_MAX: Maximum separation of charged lepton pairs, ARu. Default is 50. 
Photon specific cuts: 

• Y_G_MAX: Maximum pseudorapidity for photons. Default is 1.5. 

• PT_G_MIN : Minimum transverse momentum for photons. Default is 20 GeV. 

• RGG_MIN: Minimum separation of photon pairs, AR^^. Default is 0. 

• RGG_MAX: Maximum separation of photon pairs, AR^^. Default is 50. 
Additional cuts: 

• RJL_MIM: Minimum separation of an identified jet and a charged lepton, ARji. 
Default is 0. 

• RJG_MIM: Minimum separation of an identified jet and a photon, ARj^. Default 
is 0. 

• RLG_MIN: Minimum separation of a charged lepton and a photon, Ai?^^. Default 



• ETAJJ_MIN: Minimum required pseudorapidity gap, Arjjj, between two tagging jets 
(the two leading jets in a Pt ordering). Default is 0. 

• YSIGN : If set to . true . , the two tagging jets are required to be found in the opposite 
detector hemispheres. Default is .false.. 

• LRAPIDGAP : If set to .true . all charged leptons are required to lie between the two 
tagging jets in rapidity. Default is .false.. 

• DELY_ JL : Minimum rapidity distance of the charged leptons from the tagging jets, 
if LRAPIDGAP is set to .true.. Default is 0. 

^^These apply only to vbf and ggf processes. 



is 0. 




14 



• GRAPIDGAP : If set to .true . all photons are required to lie between the two tagging 
jets in rapidity. Default is .false.. 

• DELY_JG: Minimum rapidity distance of photons from tagging jets, if GRAPIDGAP is 
set to .true.. Default is 0. 

• MDIJ_MIM: Minimum dijet invariant mass of two tagging jets. Default is GeV. 

• MDIJ_MAX: Maximum dijet invariant mass of two tagging jets. Default is 14000 GeV. 

• JVETO: If set to .true., a central jet veto is applied. Default is .false.. 

• DELY_ JVETO : Minimum rapidity separation of a central jet from two tagging jets. 
Default is 0. 

• YMAX_VETO: Maximum pseudorapidity of a central jet. Default is 4.5. 

• PTMIN_VETO: Minimum transverse momentum of a central jet. Default is 10 GeV. 

4.7 Parameters for anomalous couplings 

Vbfnlo supports anomalous HVV couplings, where V = W,Z,'j, in both the produc- 
tion and the decay of a Higgs boson in VBF type reactions, i.e. for ProcIds 100-107. 
The anomalous HVV couplings can be parameterized in the anom_HVV.dat input file. 
Moreover, the triple and quartic anomalous gauge boson couplings for the VBF process 
pp — ^ W~^W~jj{j) are included [23]. These can be set in the anoiii_WW.dat file. 

4.7.1 anom_HVV.dat — anomalous HVV couplings 

Among the anomalous coupling input parameters, the user can choose between three 
different parameterizations. 

1. A parameterization in terms of couplings in the effective Lagrangian approach. In 
Vbfnlo effective dimension five operators are included corresponding to 

HVV HVV ^ 

^tf = ^-^HV^^V"^^ + ^i/y+y-^'^ , (1) 

where the subscript e or o refers to the CP-even or CP-odd nature of the individual 
operators [24]. 

• PARAMETRl : Parameter which switches on the effective Lagrangian parameter- 
ization Eq. dl]). The default value is .false.. 

• LAMBDAS : Mass scales Ase and A50 in units of GeV with 480 GeV chosen as 
default. 

• G5E_HWW, G5E_HZZ, G5E_HGG, G5E_HGZ : Parameters which determine the cou- 
plings g^^^ of the CP-even dimension five operators. Their default values are 
set to 0. 

• G50_HWW, G50_HZZ, G50_HGG, G50_HGZ : Parameters which determine the cou- 
plings gl^^^ of the CP-odd dimension five operators. Their default values are 
set to 0. 
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The parameterization of the anomalous couplings by the L3 Collaboration as given 
in Ref. [25]. The parameters are d, ds, Agf and Ak^. For the CP-odd operators 
only three parameters are needed. These are defined in analogy to the CP-even 
ones and can be related to the coefficients fi of the operators Cj in the effective 
Lagrangian as described in Ref. [26] in the following way: 

2 2 
"~ ^2 JWW ■> ^2 JWW ^ 



_ ?,ir?eyj ~ _ sin^^^ 



A _ 1 _ I f , f \ ~ _ '^W r 

l\K^ - K^ - I - [Jb + JW) , l^^i - Jb ' 



BB ' 

(2) 



_ z _ ^ _ !^ fw_ 
-9i i - ^2 2 

• PARAMETR2 : Parameter which switches on the above mentioned parameteriza- 
tion. The default value is .false.. 

• D.EVEW, DB.EVEN, DGIZ.EVEN, DKGAM.EVEN: Parameters which are the CP- 
even couplings in this parameterization with default 0. 

• D_ODD, DB_ODD, KGAM_ODD: Parameters which are the CP-odd couplings in 
this parameterization with default values equal to 0. 

• HVVl : Parameter which determines which anomalous HVV couplings are used 
for the run. For HVVl = 0, only the HZ'y couplings, for HVVl = 1, only the H'yy 
coupling, for HVVl = 2, only the HZZ coupling and for HVVl = 3, only the 
HWW coupling is used. If HVVl is set to 4, all possible anomalous couplings 
are used. This is also the default value. 

The parameterization of the anomalous couplings in terms of coefficients /j/A^ of 
the operators Owwt Obb, Ow and Ob and their corresponding CP-odd operators 
according to Refs. [27,28]. 

• PARAMETR3 : Parameter which switches on the parameterization stated above. 
The default value is . false . . 

• FWW.EVEM, FBB.EVEM, FW.EVEM, FB.EVEM: Parameters which represent the 
coefficients of the CP-even operators with default values equal to 0. 

• FWW_ODD , FBB_ODD , FB_ODD : Parameters which are the coefficients of the CP- 
odd operators with default values 0. 

• HVV2 : Parameter which allows to choose which anomalous HVV couplings are 
used. For HVV2 = 0, only the HZ'j coupling, for HVV2 = 1, only the H'yy 
coupling, for HVV2 = 2, only the HZZ coupling and for HVV2 = 3, only the 
HWW coupling is used. If set to 4 all possible anomalous couplings are used. 
The default value is 4. 

• TREEFAC: Parameter that multiplies the HVV tensor present in the SM La- 
grangian. Default is 1. 

• LOOPFAC: Parameter that multiplies the HZ'j and Hjj vertices induced by 
SM loops. The default is chosen to be 1. 
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Moreover, for all parameterizations two different form factors can be chosen as described 
in Refs. [24,26]. They model effective, momentum dependent HVV vertices, motivated 
from new physics entering with a large scale A at the loop level. 

A^ 

F2 = -2A'Co{qlqliqi + q2)\A'). (4) 

Here the are the momenta of the vector bosons and Cq is the scalar one-loop three 
point function in the notation of Ref. [29]. 

• FORMFACTOR: Parameter which switches on the above form factor parameterization. 
The default value is set to .false.. 

• MASS_SCALE : Characteristic mass scale of new physics A in units of GeV. The default 
value is 200 GeV. 

• FFAC : Parameter which is used to select one particular form factor out of Eqs. ((3l) 
and If FFAC = 1, the form factor Fi is used for the parameterization. FFAC = 
2 selects F2, which is also the default value. 

Finnaly, the two parameters can be used to rescale the SM HVV couplings. 

• TREEFAC : Parameter that multiplies the HVV tensor present in the SM Lagrangian. 
Default is 1. 

• LOOPFAC : Parameter that multiplies the HZ'y and Hj'j vertices induced by SM 
loops. The default is chosen to be 1. 

4.7.2 anom_WW.dat — anomalous triple and quartic gauge boson couplings 

The triple and quartic anomalous gauge boson couplings can be set in anom_WW.dat. The 
input values are the coefficients /j/A^ of the CP-even operators in the effective Lagrangian 
Obw, Odw, Owww, C>wwi Obbi Ow and Ob and their corresponding CP-odd operators 
as described in Refs. [28,30]. 

• FBW, FDW, FWWW, FWW, FBB, FW, FB: Parameters which give the values of the 
coefficients of the CP-even operators. The default values for these parameters are 0. 

• FWWt, FBWt, FBBt, FWt , FBt , FWWWt , FDWt : Parameters which are the coeffi- 
cients of the CP-odd operators. The default values for these parameters are 0. 

• OVS : Parameter which switches on the overall factor scheme. The default value is 
.false . . 

• FORMFAC : Parameter which allows to include the form factor 

1 + , (5) 

in the calculation, where A is a characteristic new physics mass scale. The default 
value is .false . . 

LAMBDA : Parameter that gives the above scale A in units of GeV. Default is 2500 GeV. 
EXPFAC: Parameter to set the exponent n in Eq. ([5]) with default value 2. 
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4.8 kk_input.dat — parameters for the Warped Higgsless model 

Vbfnlo is capable of calculating the weak boson fusion processes VV+2 jets in the 
Warped Higgsless scenario [31] at LO and NLO QCD level (see, for example, Ref. [32] for 
a phenomenological application). The model parameters can be generated by Vbfnlo 
via the input file kk_input.dat for a hard wired choice of the relevant five dimensional 
gauge parameters. The input values are 

• Switch that determines whether Vbfnlo should calculate the parameters needed 
by the model. Default is .true.. 

• Location of the UV brane for the generation of the model parameters in the Warped 
Higgsless Model. Default is R = 9.75 x 10~^, which amounts to the Kaluza-Klein 
excitations having masses of = 700 GeV, = 695 GeV, and m^/ = 718 GeV. 
Smaller values of R result in a heavier Kaluza-Klein spectrum. 

• The maximum number of Kaluza-Klein states to be included on top of the 
Standard Model bosons, which correspond to VF^^i. All states A; > 3 are 
phenomenologically irrelevant. Default is 1. 

• The maximum number of Kaluza-Klein Zk states to be included on top of the 
Standard Model Z boson, which corresponds to Zk=i. All states A; > 3 are phe- 
nomenologically irrelevant. Default is 1. 

• The maximum number of Kaluza-Klein Z^. bosons that are excitations of the Stan- 
dard Model photon Z^^q. States k > 2 are phenomenologically irrelevant. Default 
is 1. 

The explicit breaking of higher dimensional gauge invariance is balanced according to the 
description of [33,34], where also more details on the model and its implementation can 
be found. 

Vbfnlo generates the text file kk_coupl_inp.dat, which documents the calculated 
model parameters, i.e. Kaluza-Klein gauge boson masses, couplings and widths of the 
specified input parameters. This file can also be used as input file for advanced users 
who want to run the code with their own set of parameters. To that end, select .false . 
in the file kk_input.dat. Vbfnlo will then calculate the gauge boson widths on the 
basis of these parameters by the decay to the lower lying states. Information on the 
widths and on the sum rules relating the various gauge boson couplings [35] are written 
to kk_check.dat. 

4.9 ggflo.dat — general parameters for gluon fusion processes 

In Vbfnlo, the double real-emission corrections to gg — > 0, which lead to 0+ 2 jet events 
at order a^, are included. Here, can be a scalar (h, H) or pseudo-scalar (A) Higgs boson 
as in a generic two-Higgs -doublet model (2HDM) of type II. Contributions contain top- 
and bottom-quark triangles, boxes and pentagon diagrams, i.e. the full mass dependence 
of the loop induced production. Interference effects between loops with bottom and top 
quarks as well as between CP-even and CP-odd couplings of the heavy quarks are fully 
taken into account. An option to use the large top mass approximation, which works well 
for intermediate Higgs boson masses, provided that the transverse momenta of the final 
state partons are smaller than the top quark mass, is also implemented. 
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Model Description 



Toy model with settings provided by the user 



in the large top quark mass limit (nit oo) 
in the large top quark mass limit {rrit oo) 
with full mass dependence of the top quark loop 
with full mass dependence of the bottom quark loop 
with full mass dependence of the top and bottom quark 



CP-even Higgs {Hsm 
CP-odd Higgs (Hsm) 
CP-even Higgs {Hsu 
CP-even Higgs {Hsm 
CP-even Higgs {Hsm 
loop 

CP-odd Higgs (A) with full mass dependence of the top quark loop 
CP-odd Higgs (A) with full mass dependence of the bottom quark loop 
CP-odd Higgs (A) with full mass dependence of the top and bottom quark 
loop 

CP-even Higgs (h) with full mass dependence of the top quark loop 
CP-even Higgs (h) with full mass dependence of the bottom quark loop 
CP-even Higgs (h) with full mass dependence of the top and bottom quark 
loop 

CP-even Higgs (H) with full mass dependence of the top quark loop 
CP-even Higgs (H) with full mass dependence of the bottom quark loop 
CP-even Higgs (H) with full mass dependence of the top and bottom quark 
loop 



Table 10: Different models for the production of a Higgs boson plus two jets via gluon 
fusion. 



Higgs boson plus two jets production via gluon fusion requires usage of the ggf lo 
executable program rather than the vbf nlo one. Moreover, the vbf nlo . dat input file 
does not need to be modified in order to set the correct value of process ID. In order 
to be more transparent to the user, this information is automatically assumed when 
running ggf lo. However, the following additional parameters have to be adjusted in the 
ggf lo .dat file: 

• PROCESSGGF: Model ID. A summary is given in Table [lOl Default is 3. 

• SUBPRQQ : Switch for the subprocesses with quark-quark initial state. Default is set 
to .true.. 

• SUBPRQG : Switch for the subprocesses with quark-gluon initial state. Default is set 
to .true.. 

• SUBPRGG : Switch for the subprocesses with gluon-gluon initial state. Default is set 
to .true.. 

• TAN_BETA: Ratio of vacuum expectation values, tan/3 = Vu/vd where f ^ = f ^ + fj. 
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• ALPHA: Neutral Higgs boson mixing angle a, which arises when the CP-even Higgs 
boson mass matrix is diagonalized to obtain the physical CP-even Higgs boson 
states, h and H . 

5 Checks 

Extensive checks for the LO and the real emission amplitudes as well as for the total 
LO cross sections have been performed for all processes implemented in Vbfnlo. Born 
amplitudes and real emission diagrams have been compared with the fully automatically 
generated results provided by MadGraph [36]. Complete agreement has been found in 
each case. Moreover, total LO cross sections with a minimal set of cuts agree with the 
respective results obtained by MadEventEI [37,38] and HELAC-PHEGAqll [39-41], a 
completely automatic parton level event generator based on Dyson-Schwinger recursive 
equations. 

All LHA event files for the LO processes have been tested with HERWIG++0 [42], 
a general purpose Monte Carlo event generator for the simulation of hard lepton-lepton 
and hadron-hadron collisions. 

As a final and very important test, comparisons with already published results have 
been made. In Ref. [43], a tuned comparison of LO and NLO QCD results for Higgs 
boson production via vector boson fusion at the LHC has been performed. Three different 
calculations have been cross checked: Vbfnlo, the results of Refs. [44,45], and the VV2H 
program^. For the dominant t- and w-channel contributions which are implemented in 
Vbfnlo, good agreement has been found. For the triboson processes a comparison for 
the production of on-shell gauge bosons without leptonic decays has been performed with 
the results presented in Ref. [46]. Again, good agreement has been found. Results for the 
CP-odd and CP-even Higgs boson production via gluon fusion have been tested against 
FeynArt^^I [47,48]. 

6 Summary & Outlook 

Vbfnlo is a fully fiexible partonic Monte Carlo program for vector boson fusion, double 
and triple vector boson production processes at NLO QCD accuracy. The simulation 
of CP-even and CP-odd Higgs boson production in gluon fusion, associated with two 
additional jets, is implemented at leading order (for this process only, the LO starts at 
one-loop level). 

Future improvements are directed along two main lines of development: Further pro- 
cesses at NLO QCD accuracy will be included (e.g., pp — ^ WW'j and pp Wjj) and 
new features will be added to the already existing processes, such as anomalous triple 
and quartic gauge boson couplings and Kaluza-Klein excitations. Higgs production via 
gluon fusion within the generic two-Higgs-doublet model will be extended to a complete 
simulation within the Minimal Supersymmetric Standard Model (MSSM), including full 

http : / /madgraph . hep . uiuc . e^7| 

http : / /helac-phegas . web . cer n. ch/helac-phegas7[ 
http : / / pro j ects . hepf orge . org/herwig/ 
^' http : //people .web .psi . ch/ spira/vv2h/ 
http : //www . f eynarts . de/ 



20 



dependence on scalar top and bottom quarks masses. Matching the NLO QCD processes 
to a parton shower at next-to-leading logarithmic accuracy is currently in progress. 
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